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Introduction 



Within a few years another new window on the universe is expected to open 
up, with the first direct detection of gravitational waves. There is keen 
interest in observing gravitational waves directly, in order to test Einstein's 
theory of general relativity and to observe some of the most exotic objects 
in nature, like black holes. In addition, the power of gravitational wave 
observations to produce more surprises is very high. 

The gravitational wave spectrum is completely distinct and, at the same 
time, complementary to the electromagnetic spectrum. The primary emit- 
ters of electromagnetic radiation are charged elementary particles, mainly 
electrons; because of overall charge neutrality, electromagnetic radiation is 
typically emitted in small regions, with short wavelengths, and conveys direct 
information about the physical conditions of small portions of the astronom- 
ical sources. By contrast, gravitational waves are emitted by the cumulative 
mass and momentum of entire systems, so they have long wavelengths and 
convey direct information about large-scale regions. Moreover, electromag- 
netic waves couple strongly to charges and so it is easy to detect them, 
but they are also easily scattered or absorbed by material between us and 
the source; gravitational waves instead couple extremely weakly to matter, 
making them very hard to detect, but also allowing them to travel to us sub- 
stantially unaffected by intervening matter, even from the earliest moments 
of the Big Bang. 

These contrasts, and the history of serendipitous discovery in astronomy, 
suggest that electromagnetic observations may be poor predictors of the phe- 
nomena that gravitational wave detectors will eventually discover. Given that 
96% of the mass-energy of the universe carries no charge, gravitational waves 
provide a first opportunity to observe directly a major part of the universe. 
It might turn out to be as complex and interesting as the charged minor 
component, the part that we call "normal" matter. 
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One of the longstanding problems of modern gravitational physics is the 
detection of gravitational waves, for which the standard theoretical analysis 
relies upon the split of the space-time metric Qab into a background metric 
plus perturbations. However, the background need not be Minkowskian in 
several cases of physical interest, nor it has to be always a solution of the 
vacuum Einstein equations. As a consequence, we are therefore aiming to 
investigate in more detail what happens if the background space-time has a 
non-vanishing Riemann curvature. 

This issue has to be seriously considered from an experimental point of 
view since the gravitational wave detectors of new generation are also de- 
signed to investigate strong field regimes: this means that the physical situa- 
tions, where only the standard Minkowski background is taken into account, 
could be misleading to achieve self-consistent results. 

In particular, several ground-based laser interferometers have been built 
in the United States (LIGO) [DEI, Europe (VIRGO and GEO) [Sill] and 
Japan (TAMA) [5] and they are now in the data taking phase for frequency 
ranges about 10~^ kHz. However, new advanced optical configurations allow 
to reach sensitivities slightly above and below the standard quantum limit for 
free test-particles, hence we are now approaching the epoch of second [S] and 
third [7| generation of gravitational wave detectors. This fact, in principle, 
allows to investigate wide ranges of frequencies where strong field regimes or 
alternative theories of gravity can be tested [8l [9l [10] . 

Besides, the laser interferometer space antenna (LISA) [11] (which is 
mainly devoted to work in the range 10~^ ~ 10~^ Hz) should fly within 
the next decade to investigate the stochastic background of gravitational 
waves. At much lower frequencies (10~^^ Hz), cosmic microwave background 
(CMB) probes, like the forthcoming PLANCK satellite, are designed to de- 
tect also gravitational waves by measuring the CMB polarization [12] while 
millisecond pulsar timing can set interesting upper limits in the frequency 
range between 10~^ ~ 10~* Hz [T3l[Il]. At these frequencies, the large num- 
ber of millisecond pulsars detectable by the square kilometer array would 
provide a natural ensemble of clocks which can be used as multiple arms of 
a gravitational wave detector fT5\. 

This forthcoming experimental situation is intriguing, but deserves a se- 
rious theoretical analysis which cannot leave aside the rigorous investigation 
of strong field regimes and the possibility that further polarization states 
of gravitational waves could come out in such regimes. For example, if one 
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takes into account scalar-tensor theories of gravity [S] or higher-order the- 
ories [9], scalar-massive gravitons should be considered. This implies that 
the standard approach where gravitational waves are assumed as small per- 
turbations (coming only from Einstein's general relativity) on a Minkowski 
background could be totally insufficient. On the other hand, the existence 
of these further polarization modes could be a straightforward solution of 
the dark matter problem since massive gravitons could be testable cold dark 
matter candidates as discussed in [16l [T7] . 

We want to face the issue of the rigorous formulation of gravitational 
wave problem in curved backgrounds. In particular, we want to perform a 
analysis of gravitational waves in the de Sitter space-time. Achieving solu- 
tions in this maximally symmetric background could constitute the paradigm 
to investigate any curved space-time by the same techniques and could have 
interesting cosmological applications if a conformal analysis is undertaken 
as, for example, in [10], where it is shown how the amplitude of cosmological 
gravitational waves strictly depends on the cosmological background. 

Some important progress in the astronomical observations of the last ten 
years [HI |19] have led in a progressively convincing way to the surprising 
conclusion that the recent universe is dominated by an almost spatially ho- 
mogeneous exotic form of energy density to which there corresponds an effec- 
tive negative pressure. Such negative pressure acts repulsively at large scales, 
opposing itself to the gravitational attraction. It has become customary to 
characterize such energy density by the term "dark". 

The simplest and best known candidate for the "dark energy" is the cos- 
mological constant. As of today, the ACDM (Cold Dark Matter) model, 
which is obtained by adding a cosmological constant to the standard model, 
is the one which is in better agreement with the cosmological observations, 
the latter being progressively more precise. Recent data show that dark 
energy behaves as a cosmological constant within a few percent error. In 
addition, if the description provided by the ACDM model is correct, Fried- 
mann's equation shows that the remaining energy components must in the 
future progressively thin out and eventually vanish thus letting the cosmo- 
logical constant term alone survive. 

In the above scenario the de Sitter geometry, which is the homogeneous 
and isotropic solution of the vacuum Einstein equations with cosmological 
term, appears to take the double role of reference geometry of the universe, 
namely the geometry of space-time deprived of its matter and radiation con- 
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tent and of geometry that the universe approaches asymptotically. 

It is now well-known that the problem of solving vector and tensor wave 
equations in curved spacetime, motivated by physical problems such as those 
occurring in gravitational wave theory and relativistic astrophysics, is in 
general a challenge even for the modern computational resources. Within this 
framework, a striking problem is the coupled nature of the set of hyperbolic 
equations one arrives at. 

The Maxwell equations for the electromagnetic potential, supplemented 
by the Lorenz gauge condition, are decoupled and solved exactly in de Sitter 
spacetime studied in static spherical coordinates. There is no source besides 
the background. One component of the vector field is expressed, in its radial 
part, through the solution of a fourth-order ordinary differential equation 
obeying given initial conditions. The other components of the vector field are 
then found by acting with lower-order differential operators on the solution of 
the fourth-order equation (while the transverse part is decoupled and solved 
exactly from the beginning). The whole four- vector potential is eventually 
expressed through hypergeometric functions and spherical harmonics. Its 
radial part is plotted for given choices of initial conditions. 

We have thus completely succeeded in solving the homogeneous vector 
wave equation for Maxwell theory in the Lorenz gauge when a de Sitter 
spacetime is considered. The decoupling technique, analytic formulae and 
plots are completely original [20] . 

Thus, we have extended this method to the wave equation of metric 
perturbations on a de Sitter background. It is possible to show that, in 
a covariant formulation, the supplementary condition for gravitational waves 
can be described by a functional $a acting on the space of symmetric rank- 
two tensors hab (metric perturbations). For any choice of one gets a 
different realization of the invertible operator P^b'^'^ (Lichnerowicz operator) 
on metric perturbations. The basic equations of the theory read therefore as 



where the Lichnerowicz operator P^^f, results from the expansion of the 
Einstein-Hilbert action to quadratic order in the metric perturbations, sub- 
ject to $a(^) = 0. Eventually, a numerical analysis of solutions has been 
performed. 




0, 



$a(/i) 
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However, we want to solve explicitly the Einstein equations for metric 
perturbations on a de Sitter background. Thus, we consider the vacuum 
Einstein equations with cosmological constant A, i.e. 

Rab - ^Rgab + ^9 ah = 0. 

If we introduce gab = lab + e^afe, where e is a dimensionless parame- 
ter which controls the perturbation, we get a coupled system of differential 
equations to first-order in the metric perturbation hah- At this stage, using 
the Regge- Wheeler gauge, we can solve this system exactly in terms of the 
Heun general functions [2T] . 
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Chapter 1 



Gravitational waves in 
de Sitter space-time 

The non-linearity of the gravitational field in general relativity is one of 
its most characteristic properties and it is likely that at least some of the 
crucial properties of the field show themselves only through the non-linear 
terms. Moreover, it is never entirely clear whether solutions derived by the 
usual method of linear approximation necessarily correspond in every case 
to exact solutions. 

General relativity is a peculiarly complete theory and may not give sen- 
sible solutions for situations too far removed from what is physically reason- 
able. The simplest field due to a finite source is spherically symmetrical but 
Birkhoff's theorem shows that a spherically symmetrical empty-space field is 
necessarily static. 

Therefore there cannot be truly spherically symmetrical waves and thus 
any description of radiation from a finite system must necessarily involve 
three coordinates significantly. This enormously complicates the mathemat- 
ical difficulties and thus one has to make use of methods of approximation. 

The standard theoretical analysis relies upon the split of the space-time 
metric gab into background plus perturbations, that is (unlike the contents 
section, here we don't write explicitily the e parameter) 

gab — lab + hab, 

where '-fab is the background Lorentzian metric, often taken to be of the 
Minkowski form riab, while the symmetric tensor field hab describes perturba- 
tions about Jab- The background jab needs not to be Minkowskian in several 
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cases of physical interest, nor it has to be always a solution of the vacuum 
Einstein equations. As a consequence, we are therefore aiming to investigate 
in more detail what happens if the background space-time (M, ^yab) has a 
non-vanishing Riemann curvature. 

In this work, we perform a analysis of gravitational waves in de Sitter 
space-time. 

1.1 Einstein's equations and de Sitter space- 
time 

Any space-time metric satisfies Einstein's field equations 

Rab - ^Rgab + ^9ab = ^TlTab, (1.1) 

where A is the cosmological constant. We shall use c — 1 and units of mass in 
which G ~ 1 (geometric units). Since both sides are symmetric, these form 
a set of ten coupled non-linear partial differential equations in the metric 
tensor components and its first and second-order derivatives. However, due to 
the so-called Bianchi identity, the covariant divergence of each side vanishes 
identically, that is, 

Vb ^i?"^ - ^Rg"'' + Ag"^^ = (1.2) 

and 

VbT"'' = 0, (1.3) 

hold independent of the field equations. Thus the field equations really pro- 
vide only six independent differential equations for the metric. This is in fact 
the correct number of equations needed to determine the space-time, since 
four of ten components of the metric can be given arbitrary values by use of 
the four degrees of freedom associated with a coordinate transformation. 

The space-time metrics of constant curvature are locally characterized by 
the condition ^ 

Rabcd — J^Rigacgbd — gadgbc)- (1-4) 

and this equation is equivalent to 

Rab - \Rgab = 0, (1.5) 
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Thus, the Riemann tensor is determined by the Ricci scalar it! alone and 
the Einstein tensor becomes 

Rab — -^Rgab = —-^RQab- (1-6) 

One can therefore regard these spaces as solutions of the field equations 
for an empty space with A = The space of constant curvature with 
R — Q is Minkowski space-time. The space for > is de Sitter space-time, 
which has the topology R^ x S^. It is easiest visualized as the hyperboloid 
in five-dimensional Minkowski space given by 

- + {x'f + {xY + {x'f + {xy = p (1.7) 
where A is related to Bubble's constant, Hq, by 

H'o - f ■ (1-8) 
In the standard spherical coordinates, {t, r, 9, 0), one has 



x^ = r sin 9 cos 0, 
x"^ — r sin 9 sin 0, 
x^ — rcos9, 



= - ^^cosh(m). 



^° = W;^ -^^sinh(m). (1.9) 

Thus, the metric becomes 

ds^ = -fdt^ + jdr^ + r\d9'^ + sin^ ^^0^), (1.10) 

where 

/ = 1 - ny. (1.11) 

Now, consider de Sitter metric in Gaussian confomally flat umbilical co- 
ordinates 

ds^ = -dt^ + e2^°*[dx2 + dy^ + dz^] , (1.12) 
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This metric satisfies matter-free Einstein's equations witli a non-vanisliing 
cosmological constant A sucli tliat Hq = A/3. Moreover, tlie timelike unit 
normal vector field to the t =constant hypersurfaces 

n = dt (1.13) 

form a geodesic and irrotational congruence; the 3-metric induced on the 
t =constant hypersurfaces resuhs conformally flat: 

9ab = e^^'^'Sab, (a, 6 =1,2,3); (1.14) 

finally the extrinsic curvature of these hypersurfaces is 

K{n)ab = -Hogab • (1-15) 

An orthonormal frame associated with n is given by 

n = dt, ea = e-^«U, (a = 1,2, 3) (1.16) 

For this metric the geodesic equations can be integrated exactly. In fact, 
they reduce to: 

^ = C^e-2^°* , (^X = -e + C^e-^^"* , (1.17) 



dA ' \dX J 

where the parameter e = 0,1,-1 for null, spacelike (with proper length 
parametrization, say X = s) and timelike (with proper time parametrization, 
say A = r) geodesies respectively, and C\ i = 1,2,3 are constants with 

The physical components of the tangent vector to the geodesies with 
respect to the frame fll.161) result then in 



n H e- 

V-e + C^e-'^Hot * 



:i.l8) 



It is convenient to discuss the three cases e = —1, 0, 1 separately, denoting 
the three different tangent vectors by f/{-i) = U ^ U{q) = P and f/(i) = T, 
respectively. For timelike geodesies we have 



U = -f{U,n)[n + u{U,n)i>{U,n)] 



cosh a(t)n + sinh a(t)— ej. 



cosha(t) = Vl + C2e-2^^o*, (1.19) 
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identifying the speed 



u^U, n) = tanha;(t) 



Ce 



-Hot 



as well as its direction (unit spacelike vector) 



i){U,n) 
For null geodesies we have 

P = E{P, n) [n + i>{P, n)] = Ce"^"* 



C 



;i.2o) 



identifying the relative energy 

E(P, n) = Ce-^o* . 



z>(P, nf 



9i 
'C 



:i.2ii 



1.22) 



;i.23) 



1.1.1 Null geodesies 



Let us consider first the null case. The general solution of Eq. (I1.17P is given 
by 

e« = i/„CA + c,. .. = -|l,^^^ + 4. (1,24) 

The integration constants Ci, can be chosen in such a way that x"(A = 
0) = Xq, whence 



-Hoto 



HoC 



;i.25) 



so that the solution fll.24p becomes 



t = \n[HoCX + e 
-no 



Hotoi 



X 



1 



HnCX + e^o*o 



The latter equation can also be cast in the form 



being 



x' = xi- —R{t,to) 



-no 



+ Xq . 

(1.26) 
(1.27) 
(1.28) 
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1.1.2 Timelike geodesies 



Let us consider now the timelike case. The general solution of Eq. (11 .1711 is 
given by 



e^«* = Csinh(i7or-ci) , x' 



coth{HoT - ci) + ci . (1.29) 



The integration constants ci, can be chosen in such a way that a;"(A 
0) = Xq, whence 



Ci = — arcsinh 



C 



1 ^2 — Xg + 



1 + C^e-^^oto 



1.30) 



so that the solution (I1.29P becomes 



X 



C sinh I HnT + arcsinh 



Vl + C2e-2^^o*o 



coth ( HqT + arcsinh 



+ Xr,. 



IT 



;i.3i) 



The latter equation can also be written as 
C 



X 



- 7F7^ fv^l + C'e~2^ot - v^l + C^e-^^'^A . (1.32) 



1.2 Conformal form of de Sitter metric 

It is well known that the de Sitter metric can be written as conformal to the 
Minkowski metric 

-2 



l + ^{x' + y' + z'-t'] 



-dt^ + dx^ + dy^ + dz"^) . (1.33) 



However, the explicit coordinate transformation allowing to cast the met- 
ric ( I1.12P in the previous form is somehow hidden in the literature. First of 
all introduce standard polar coordinates 



x = p sin 6* cos (p , y = p sin 9 sin (p , z = p cos 9 . 



1.34) 
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The line element fll.l2p thus takes the form 



;i.35) 



Applying then the following coordinate transformation 



t = — In 
-no 



P 



Re-^o^ 



e = e 



gives 



;i.36) 



:i - HlR^)dr^ + ^^1^ + R^{d9' + sin^ --^^ 



:i.37) 



The further transformation 



In 



2Ho 
finally gets 

ds' = 



Hlp^-{H,i-2f 



R 



1 + f (p2_p) 



9 = 6 



;i.38) 



[-dr + dp' + p\de' + sin' ed(p')] , (1.39) 



which reduces to the line element (11.331) once the Cartesian coordinates are 
restored by using standard relations as in Eq. (I1.34p . 

By combining the transformations (I1.36P and (I1.38P we get 



In 



P 
6 



Hlp^-{H,i-2f 
Hlp^-{H,i+2Y 

p Hlp^-{H,t-2Y 
/7oV- W+2)2 



^oP 



1 + f (p2_p)]2^ 



1 + f (p2_p) 



2^2 



H-^P 



[l + i^(p2_p)p 



-1/2 



;i.4o) 



which allows to pass directly from the metric (I1.35P to the conformal one 
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1.3 Preservation of the de Donder supple- 
mentary condition 

In classical gauge theory with space-time metric of Lorentzian signature, the 
gauge-fixing (also called "supplementary" condition) leads to a convenient 
form of the field equation for thepotential. For example, for classical elec- 
trodynamics in the Lorenz gaugqj, the wave equation reduces to equation 
(see Chapter 2) 

\jAb-R:A, = Q. 
However, while Maxwell Lagrangian 

Lem = -\F,bF''\ (1.41) 
is invariant under gauge transformations 

Al = A, + V J, (1.42) 
where / is a freely specifiable function of class C^, the Lorenz gauge 

$(A) = V% = 0, 

as well as any other admissible gauge, is not invariant under (ll.42p . Never- 
theless, to achieve the desired wave equation on Af,, it is rather important 
to make sure that both Ai, and the gauge-transformed potential AI obey the 
same gauge-fixing condition, i.e. [22] 

^{A) = 0, $(^0 = 0. (1.43) 

A more general situation, here not considered, is instead the case when 
only the gauge-transformed potential obeys the gauge-fixing condition, i.e. 

m 

^{A) ^ 0, $(A^) = 0. (1.44) 

The counterpart of fll.43p for pure gravity is the well known problem of 
imposing a gauge on metric perturbations and then requiring its invariance 

-'^In [24|, the author L. Lorenz, who was studying the identity of the vibrations of Ught 
with electrical currents, built a set of retarded potentials for electrodynamics which, with 
hindsight, can be said to satisfy the gauge condition W^Ai, = 0, which therefore should 
not be ascribed to H. Lorentz. 
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under infinitesimal diffeomorphisms. It is straightforward to sliow that, in a 
covariant formulation, the supplementary condition for gravitational waves 
can be described by a functional acting on the space of symmetric rank- 
two tensors hat- For any choice of ^a, one gets a different realization of the 
invertible operator P^f,'^'^ on metric perturbations. The basic equations of the 
theory read therefore as 

Pab^'^^cd = 0, 

Mh) = 0, 

where P^b'^'^ results from the expansion of the action functional to quadratic 
order in the metric perturbations. In general relativity, if one wants to obtain 
the standard covariant wave operator on metric perturbations, this is taken 
to be of the de Bonder type [25] 

Mh) = V (^Kt - ^labh^ , (1.45) 

where h = 'y'^'^hcd and denotes covariant derivative with respect to the 
background metric '^ab- Under infinitesimal space-time diffeomorphisms, the 
metric perturbations suffer the variation (the round brackets denoting sym- 
metrization) 

Shab = V(a iPb), (1-46) 

where ipb is a covector, with associated one-form iphdx^ and vector field '■p"'-^ 
(having set ip"" = 'y^'^fb, which results from the isomorphism between tangent 
and cotangent space to the background space-time, that turns covectors into 
vectors, or the other way around). The change suffered from the de Bonder 
gauge in f ll.45|) when metric perturbations are varied according to f ll.46p is 
then found to be 

6^a{h) = - {6,'\J + Ra) ^t, (1.47) 

where Q is the standard d'Alembert operator in curved space-time, i.e. 

□ = (1.48) 

By virtue of Eqs. fll.45p and (11.471) . if the de Bonder gauge was originally 
satisfied, it is preserved under space-time diffeomorphisms if and only if ipb 
solves the equation 6^a{h) = 0, that is 



-\Z\^a = Ra^b- 



(1.49) 



20 CHAPTER 1 . GRAVITATIONAL WAVES INDE SITTER SPACE- TIME 



At this stage, adding ipi, to both sides of (ll.49p . one has 



Pj' V, = 2R^ 



:i.50) 



where 



Pa'^-m+R^ (1.51) 

is the standard gauge-field operator in the Lorenz gauge (see Chapter 2). 
Thus, we can solve Eq. f ll.SOp in the form 

y., = ypf + 2P,'^ (1.52) 

where (p^^ is a solution of the homogeneous wave equation [2S] 

0, (1.53) 

while PJ^ is the inverse operator of Pj', satisfying 



1.54) 



The operator P^ is an integral operator with kernel given by the mass- 
less spin-1 Green function Gab{x,x') = Gat'- The latter can be chosen, for 
example, to be of the Feynman type, i.e. that solution of the equation (see 
Appendix A for the notation) 

6{x, x') 



Pa Ghd — 9a 



'-1 



having the asymptotic expansion as cr — [271 EH] 



Gnh' 



A-^^ + Kt' log(a + \e) + 



o + \e 



ah' 



;i.55) 



1.56) 



where (j{x, x') is the Ruse-Synge world function [291 |30l [3T] , equal to half the 
square of the geodesic distance /i between the points x and x'. 



1.4 Massless Green functions in de Sitter space- 
time 

This general scheme can be completely implemented in the relevant case [32] 
of de Sitter space where, relying upon the work in [33], we know that the 
massless spin-1 Green function reads as 



Gab' = (y{fj,)gab' + /3{fj,)nanb', 



(1.57) 
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where x') = ^y2cr{x,x') is the geodesic distance between x and x', 
n'^{x,x') and n"''{x,x') are the unit tangents to the geodesic at x and x', 
respectively, for which 

na{x,x') = Vafi{x,x'), 

na'{x,x') =Va'fi{x,x'), (1.58) 
while, in terms of the new variable 



z = - \ 1 + cos 



p 

the coefficient functions a and /3 are given by [33] 



487r2p2 



3 1 

1-z) z 



2 1 , , . 
- + - log(l - 2;) 



247r2p^ 



1 - 



1 



1 



log(l - z) 



[1.59) 



;i.60) 



:i.61) 



Strictly speaking, the formulae fll.60p - fll.6ip are first derived in the Eu- 
clidean de Sitter space. In the Lorentzian de Sitter space-time M which is 
what we are interested in, one can define the set 



Jj. = {x' G M : 3 geodesic from x to x'} . (1.62) 

Moreover, it is well-known that M can be viewed as an hyperboloid 
imbedded in fiat space, i.e. as the set of points Y"- E R""^^ such that 

where rjab = diag(— 1, 1, 1, 1), so that its induced metric reads as 

ds^ = r]abdY''dY\ 
As is stressed in Ref. |33], the relation 



;i.63) 
;i.64) 



z{x, x') 



1 + 



7]abY''{x)Y\x') 



:i.65) 



is well defined both inside and outside J^, and it is an analytic function of 
the coordinates Y"'. Thus, Eq. f ll.65p makes it possible to define z{x^x') 
everywhere on de Sitter, and one can define the geodesic distance 

/i(x, x') = 2pcos^^(v^) (1.66) 

as the limiting value [33J above the standard branch cut of cos~^. Along 
similar lines, the equations defining na,na' and gab' have right-hand sides 
which are analytic functions of the coordinates F", and are hence well defined 
everywhere on Lorentzian de Sitter space-time [33]. 
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1.5 Evaluation of the kernel 



In a de Sitter background the Ricci tensor is proportional to the metric 
through the cosmological constant: Rat = ^Qab-, and hence the formulae 
( ll.52p . f ll.57p . f ll.60p and fll.6ip lead to the following explicit expression for 
the solution of the inhomogeneous wave equation fll.SOp : 



;i.67) 



+ /3(z(/i(a;, x')))«c^"' (pa'{x')^J --i{x')d'^x' , 
where, from Eq. (I1.66p . 



/x(a;, x') = 2pcos 



-1 



^ 1 + 



7]abY''{x)Yb{x') 



:i.68) 



while Eqs. fll.60p and f ll.6ip should be exploited to express a and (3, bearing 
in mind Eq. (12.60 jointly with 



z{x, x') 



1 + COS 



/i(x, x') 

p 



:i.69) 



Moreover, the bivector in the integrand fll.67p is given by 



9a 



1 

psin(/i/p)' 



;i.70) 



The right-hand side of the formula expressing gj^' is an analytic function 
of the coordinates Y°- and is therefore well defined everywhere on de Sitter 
[55] . The integral on the right-hand side of Eq. ( I1.67P can be conveniently 
expressed the form 



Ux) 



a(2)C-^(p)VcV'''/i 



+ {P{z)-a{z)){V,^^){V'''^l)\ipa'{x')^^^)d^x', (1.71) 
with a and (3 — a given by (cf. f ll.60p and f ll.6ip ) 

;i + 2z) 



a{z) 



z(l-z] 



+ — log(l-2;) 



;i.72) 
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/3{z) — a{z) 



-S + 2z- 2z') 3 



487r2p^ 

Equation fll.67p is therefore an integral equation reading as 



with unbounded kernel given by 



KJ" =2 



This kernel is indeed unbounded by virtue of the limits 



;i.73) 



(1.74) 



a(^)C-i(;u)VeV"V + - a(^))(Ve/x)(V'^V)l • (1-75) 



487r^p^ lim za{z 



2->0 



487rV \im(l-z)a(z) = 1, 

2-)-l 



487rV lim z{l3{z) - a{z)) = --, 
z— >o 2 

487rV^ lim(l - z)((3(z) - a(z)) = -3. 



(1.76) 

(1.77) 

(1.78) 
(1.79) 



At this stage, we can exploit fll.65p and fll.75p to re-express the kernel in 
the form 



{V,z){V-'z) 



+ 



247rV(l - z 
(VcV"'z 



2 + ( -3 + ^(1 + 2^: 



+ 4 log(l - z] 



z{l — z) z"^ 



z{l + 2z) 



+ — log 1-2; 



zil-z) z^ 



;i.8o) 



Note now that Lpf\x) in Eq. fll.74p . being a solution of the homogeneous 
vector wave equation (11.531) . admits the Huygens principle representation [27] 



'S' L ' J 



where 



cb' 



[1.81) 



(1.82) 



Gcb';m' — -^Kcb'-m'- 



:i.83) 



Unlike the work in |27], we here advocate the use of the Green function 
(11.571) rather than the sum, over all distinct geodesies between x and x', of 
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the Hadamard functions. To lowest order in the cosmological constant A, 
Eq. fjl.Sip may be used to approximate the desired solution of Eq. f ll.74p in 
the form 

ip,{x) = ip'-^\x) + A [ K/Vi°^ V^^tM^^V + 0(A2). (1.84) 



Omitting indices for simplicity, the general algorithm for solving Eq. 
fll.74[) . here re- written in the form 



(p = (p 



+A J Kip, (1.85) 



would be instead 



= <^(o)+A j ir<^(o), (1.86) 
<^2 = (^(°) + A J Kipi = (^(°) +aJ Kif^'^^ +A^ J J KKip^'^\ (1.87) 



+ ^A^ I ... j K^^^^\ (1.88) 

ip = lim ipn- (1.89) 

In this Chapter, we have seen that when the de Bonder gauge is imposed, 
its preservation under infinitesimal space-time diffeomorphisms is guaranteed 
if and only if the associated covector is ruled by a second-order hyperbolic 
operator which is the classical counterpart of the ghost operator in quantum 
gravity and the vector wave equation (11.491) is been studied by using an 
integral representation. 

However, a different approach is viable, that is, through a solution by 
factorization of a hyperbolic equation. In fact, in the equation fll.491) the Ricci 
term has opposite sign with respect to the wave equation for Maxwell theory, 
in the Lorenz gauge. Thus, we are interested in the following generalized wave 
equation: 

- □X, + ei?„% = 0, (1.90) 

where e = ±1. In particular, e = 1 corresponds to studying the Maxwell 
vector wave equation, whereas e = — 1 provides our consistency equation 
( 09|) . 

In the Chapter 2, by virtue of the spherical symmetry of de Sitter space- 
time, these equations should be conveniently written by using the expansion 
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of X in vector harmonics. In this way, we solve the Maxwell equation in de 
Sitter space-time [20] and, at this stage, it is possible to relate the solutions 
of the two problems [25] . 
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Chapter 2 

The vector wave equation 



It is by now well known that the problem of solving vector and tensor wave 
equations in curved space-time is in general a challenge even for the modern 
computational resources. On using the Maxwell action functional 



S^-\ I F^.F'^'^ d^x (2.1) 
4 -Im 



where Fah is the electromagnetic field tensor, one gets the wave operator P^, 
that is, 

^a'=-C□+^a+V„V^ (2.2) 

but jointly with the Lorenz gauge condition 

V% = 0, (2.3) 

we have 

^a%=Han+^a)A- (2-4) 

Thus, in vacuum, the coupled equations for the electromagnetic potential 

are 

{-8aU+K)A,^Q (2.5) 

and eventually 

□ Afc - R,'A, = 0. (2.6) 

We note that in the quantum case, one has, even on using the Lorenz 
gauge condition, the wave operator 

Pa' = -CD + Ra + - ^) ^"^^ (2.7) 

27 
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but, following Feynman, we put a = 1 and one has 



(2. 



that is, at least formally, the same wave operator of the classical case. 

A deep link exists between classical and quantum theory, since in the 
latter, the one-loop analysis depends on the functional determinant of the 
operator Pj'. 

At this point, we want to study the vector wave equation in de Sitter 
space-time with static spherical coordinates, so that the line element is fll.lOp . 
The vector field X solving the vector wave equation can be expanded in 
spherical harmonics according to 



X = Yim{9)e 



—i{ujt—m(j>) 



fo{r)dt + h{r)dr 



+ e 



sin 6'' 



de 



de 



+ ie 



-i{ujt—mij>) 



dYi ~ 
-r sin 6^/2 (r)—^ + mf^{r)Yira{e) 



(2.9) 



where Yim{e) is the ^-dependent part of the spherical harmonics Yir^ 
solution of the equation 



^2 ^ d 



m 



+ L 



sm 



YiUO) = 0, 



(2.10) 



with L = /( / + 1). As one has shown in [25], these equations lead to a system 
of coupled ordinary differential equations for the functions /o, /i, /s, besides 
a decoupled equation for /2 (/2 being related to the transverse part of X). 
The equation for /2(r) can be easily integrated in terms of hypergeometric 
functions. In fact, assuming 



/2(r) = /-^^/(2^V(r), 



(2.11) 



the resulting equation for ip reads as 

-— = 2iH r -1 + uHr — 

dr'' rj V / dr 



2 2 
u r 



L 



2i7V2 + SiuHr'^ 
(2.12) 
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Thus, the solution f2{r) is of the form 

/2(r) = r^/' [f/ir'F (^a_, ^ + /; H\'^ +U,r-^-'F (a+, a_; ^ - /; i/V^ 

(2.13) 

where 

Q = ^,a± = — {2in-3-2l±l). (2.14) 
H 4 

At this stage, however, the problem remained of solving explicitly also 
for /o(r), /i(r), fs^r) in the expansion f l2.9p . For this purpose, we derive the 
decoupling procedure for such modes in de Sitter and we write explicitly the 
decoupled equations. Eventually, we solve explicitly for fo, fi, fs in terms 
of hypergeometric functions and we plot such solutions for suitable initial 
conditions. 



2.1 Coupled modes 

Unlike /2, the functions /o, /i and /a obey instead a coupled set, given by 
Eqs. (54), (55), (57) of [25], which are here written, more conveniently, in 
matrix form as (we set e = 1 in the Eqs. of [25], which corresponds to 
studying the vector wave equation fl2.6p ) 

a 

Pi r-'f'L^ I I I = 0, (2.15) 





having defined 



a = — (2.16) 
/3 ^ -, (2.17) 



(2.18) 
(2.19) 
(2.20) 







^ d 


Po = 




+ Qi-r- 

dr 






^ d 


Pi = 


dr'^ 


dr 




d^ 




Ps = 


dr^ 


+ Q5-r 

dr 
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Qi ^ /? = -, (2.21) 
r 

- ^-^^ (2.22) 

Qs = -fl-^T^ (2-23) 



^ 5^-(-'-^)7' (2-) 

= (2.25) 



«. - 7^ -^7- (2.26) 

With our notation, the three equations resulting from f l2.15p can be writ- 
ten as 



Pofo = (2.27) 
Pi/i = -^/o-^/3, (2.28) 
Psh = -/3fi. (2.29) 



2.2 Decoupled equations 



We now express /i from Eq. fl2.27p and we insert it into Eq. fl2.28p . i.e. 



Next, we exploit the Lorenz gauge condition (12.31) . i.e. |25] 



(2.30) 



f-''^U^Pf) 2r(l-2/) / 1 \ /]gr^ 



and from Eqs. f l2.30l) and fl2.31l) we obtain, on defining the new independent 
and dimensionless variable x = rH, the following fourth-order equation for 

u 



d 



/o(x) = 0, (2.32) 
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where 



(2.33) 



- 1)4' 

k{x) = L(L-2) + 2L(2-L-^]2)x2+b^ + 4^]2 + L(L + 2(^]2-l))jx^ (2.34) 

4(^]2 + + 6x2) 



^2(3;) 



Ki[X 

2 



x{x'^ — ly ' 

L + (fi2 + L - U)x^ + 18x^ 



K3(x) 



^^(x^ — 1)2 
4(-l + 3x2) 



(2.35) 



(2.36) 



(2.37) 



x{x'^ — 1) 

Eventually, fi and F3 = iJ/a can be obtained from Eqs. f l2.27p and f l2.3ip . 



I.e. 



i {I — x^) / d 



- 



2 d 
+ — r- + 



fi2 



Fsix) 



2fi X ydx"^ X dx (1 — x2)2 x2(l — x^) 
x2(l - x2) 2x(2x2 - 1) 



/o(x), (2.38) 



L .X L j^^(^) + 4(r^^°^^)- ^2-^^) 

Our /i and /s are purely imaginary, which means we are eventually going 
to take their imaginary part only. Moreover, as a consistency check, Eqs. 
(12381) and (^Wf have been found to agree with Eq. ([229]), i.e. (Km is 
then identically satisfied. 



2.3 Exact solutions 

Equation (12.321) has four linearly independent integrals, so that its general 
solution involves four coefficients of linear combination Ci, C2, C3, C4, accord- 
ing to (hereafter, F is the hypergeometric function already used in (12.131) ) 

foix) = CiX-^-'(l-xr^-F(-^l]-i,-^l]+i-^;i-/;x^ 

+ c,x\i - x^r^^F [-Ui + -^^ + ^ + ^ + 

+ C4x'(l - x2)-t^F (-'-n + 1 + ^> + ^ + ^ + ^')2-40) 
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Regularity at the origin (recall that x = should be included, since the event 
horizon for an observer situated at x = is given by x = 1 [35J) implies that 
Ci = C2 = 0, and hence, on defining 

ai = ~n + bi = -'-n + ^+^-, d, = ^ + i, (2.41) 

we now re-express the regular solution in the form (the points x = 0, ±1 
being regular singular points of the equation 02.321) satisfied by /o) 

/o(x) = x\l-x^)-^^^\CsF{ai,bi;di;x'^)+CiF{ai + l,bi-l;dr,x'^)\, (2.42) 

where the second term on the right-hand side of 02.421) can be obtained from 
the first through the replacements 

C3 — 7- C4, ai — 7- ai + 1, 61 —> 61 — 1 

and the series expressing the two hypergeometric functions are conditionally 
convergent, because they satisfy Re{c — a — b) = iQ, with 

a = oi, oi -|- 1; 6 = 61, 61 — 1; c = di. 



Last, we exploit the identity 



^F(a, b; c; z) = —F{a + 1, 6 + 1; c + 1; z) 
dz c 



(2.43) 



to find, in the formula 02.38P for fi{x 
d 



dx 



Mx) 



+ — — x'+'(l 
di 



1(1 - x'^) + i^lx'^ F{ai, bi, di, x^) 
F{ai + 1,61 + l;di + l;x^) 



jCs C4, ai ai + 1, 61 61 - l|. 



(2.44) 



It is then straightforward, although tedious, to obtain the second deriva- 
tive of /o (see Eq. OB.ip of the Appendix B) in the equation for fi and the 
third derivative of /o in the formula 02.391) for Hf^. The results are exploited 
to plot the solutions. 

In general, for given initial conditions at x G [0, 1[, one can evaluate C3 
and C4 from 

fo{x) = y, foix) = y', 
i.e. Cs = Csiy,y'),C, = C,iy,y'). 
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Figure 2.1: Regular solution I^M) for /o with C3 = 0, C4 = 1, / = 2, . . . , 10 
with Q = 2 (left figure) and Q = A (right figure). Increasing values of / 
correspond to more peaked curves on the right part of the plots. 




Figure 2.2: Regular solution (^M) for /o with C3 = 1,0^ = 0,1 = 2, . . . ,10 
with Q = 2 (left figure) and Q = 4 (right figure). Increasing values of I 
correspond to more peaked curves on the right part of the plots. 

2.4 Plot of the solutions 

To plot the solutions, we begin with /q as given by (I2.42p . which is real- valued 
despite the many i factors occurring therein. Figures (2.1) to (2.3) describe 
the solutions for the two choices C3 = 0, C4 = 1 or the other way around and 
various values of I and Q. 

We next plot fi/i and F^/i = Hf^/i by relying upon fl2.38p and (12.391) . 
As far as we can see, all solutions blow up at the event horizon, corresponding 
to X = 1, since there are no static solutions of the wave equation which are 
regular inside and on the event horizon other than the constant one [35]. 
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0.2 0.4 O.S 0.8 



Figure 2.3: Regular solution (EH for /o with C3 = 0, C4 = 1, / = 2 and 
Q = 10. 




Figure 2.4: Regular solution fl238|) for fi/i with = 1, ^4 = (left figure) 
and C3 = 0, C4 = 1 (right figure) for / = 2 . . . , 10 and 1] = 4. 




Figure 2.5: Regular solution fl2.39p for F^/i with C3 = 1, C4 = (left figure) 
and C3 = 0, C4 = 1 (right figure) for Z = 4 and = 4. 
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2.5 Special cases: I = and / = 1 

We list here the main equations for / = 0, 1 for completeness. 

2.5.1 The case / = 0,m = 

In this case we have Y{6) = (An)'^/'^ — constant and the only surviving 
functions are fo and /i. The main equations reduce then to 

d^fo _ 2dfo -3e + 3ex^ + 3-3x\ Ml^ 

dx^~ xdx {x^-lf ^o + ^2_i/i 1^-4^) 

and the Lorenz gauge condition which now becomes 

a/„^,(.^-l)^|^ ^'(-'-f-'-^ /.. (2.46) 

These equations can be easily separated and explicitly solved in terms of 
hypergeometric functions. 

2.5.2 The case / = l,m = 0, 1 

In this case we have 

Y = \ —cos9 / = l,m = 0, 

V 47r 

Y = -W-^sin^ / = i,m=l. (2.47) 

V oTT 

However, due to the spherical symmetry of the background the equations 
for both cases I — l,m — and Z = 1, m = 1 do coincide. We have 

dVo _ _ 2 d/o _ ^'^x^ + 3ex^ + 5^2 - 3x^ - Sex"^ - 2 ^ ^ 2inx ^ 
dx'^ X dx x^{x'^ ~ 1)^ — 1 ' 

dVi _ 2(3x2 - 1) d/i 3ex^ + x^ - 3ea;2 + 3^2 + Q V - 4 2inx 

+ (x2-l)3^° 



dx'^ 


x(a;2 — 1) 
1 4F3 


dx 


x2(x2 


- 1)2 




x^{x'^ — 1) 










2(2x2 - 1) 


df2 n'^x^ - X* 


-3ex2 


- 2 + 3ex^ + 3x2 


dx'^ 


x(x2 — 1) 


dx 


x2(x2 


- 1)2 


d'Fs 

1 o 


2x dF3 

n t 1 


Q^x^ + 3ex^ H 


- _ 

o / n 


3x^ - 3ex2 - 2 ^ 



/2, 

2 

2 1 J 2 1 " -^l" (2-48) 

x^ — 1 dx x^(x^ — ly X 
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To this set one must add the Lorenz gauge condition, which now reads 

Q/o = i{x' - 1)^£ + ^^^"^^"^^ [^3 + hx{2x' - 1)]. (2.49) 

Once more, the detailed discussion of this case can be performed repeating 
exactly the same steps done in the general case. 



Chapter 3 



The tensor wave equation 



3.1 Metric perturbations in de Sitter space- 
time 

One of the longstanding problems of modern gravitational physics is the 
detection of gravitational waves, for which the standard theoretical analy- 
sis relies upon the split of the space-time metric gab into "background plus 
perturbations", i.e. 

gab = lab + hab (3.1) 

where is the background Lorentzian metric, often taken to be of the 
Minkowski form rjab-, while the symmetric tensor field hab describes pertur- 
bations about 7ab. However, the background need not be Minkowskian 
in several cases of physical interest. As a consequence, we are therefore aim- 
ing to investigate in more detail what happens if the background space-time 
(M, 7afe) has a non-vanishing Riemann curvature. In particular, we perform 
an analysis of gravitational waves on a de Sitter space-time. 

It is straightforward to show that, in a covariant formulation, the supple- 
mentary condition for gravitational waves can be described by a functional 
$a acting on the space of symmetric rank-two tensors hab occurring in Eq. 
(13. ip . For any choice of $a, one gets a different realization of the invertible 
operator Pab'^'^ (Lichnerowicz operator) on metric perturbations. The basic 
equations of the theory read therefore as (see also Chapter 1) 

Pab^'^^cd = 0, 

Mh) = 0, (3.2) 
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where the Lichnerowicz operator P^^'^'^ results from the expansion of the 
Einstein-Hilbert action to quadratic order in the metric perturbations, sub- 
ject to ^aih) = 0. 

Consider de Sitter metric in standard spherical coordinates fll.lOp . It sat- 
isfies the vacuum Einstein equations with non-vanishing cosmological con- 
stant A such that = |-. Moreover, the time-hke unit normal vector fields 
n to the t = constant hypersurfaces, 

n = dt, (3.3) 

form a geodesic and irrotational congruence. The 3-metric induced on the 
t = constant hypersurfaces turns out to be conformally fiat and the mixed 
form of the space-time Ricci tensor is simply given by 

= 3H^6a'. (3.4) 

The aim of the present chapter consists in studying the invertible wave 
operator P^f^'^'^ on metric perturbations. On considering the DeWitt super- 
metric 

j^abcd ^ ^a(c ryO-b^cd^ ^3_g-j 

the de Bonder gauge in Eq. 03.51) can be re-expressed in the form 

$„(/i) = E^'^'^Vkhai, (3.6) 

and the resulting Lichnerowicz operator on metric perturbations, obtained 
by expansion of the Einstein-Hilbert action to quadratic order in hah, subject 
to ^a{h) = 0, reads as (see [36], [37], [38] and Appendix C) 

Pab'" = Ej'i- B+R)- '^Ej^R^Hfl'^ - Ej'Rf - EJ^Rf. (3.7) 

A wave equation for metric perturbations is therefore given by 

Pah'^h.a = 0. (3.8) 

In de Sitter space-time the Lichnerowicz operator becomes 

Pah'' = Eat"' (^-D + ^R^ + |7a67 

SO that the wave equation then reads as 

= Pab^'Kd = □ + ^i?) Kb + ^labh, (3.10) 



cd 



(3.9) 
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or 



implying also 



2 \ R 

+ Kb - -^labh = 0, 



\3 + lR]h-'^Rh = -\3h = o, 



after contraction with 7 



ab 



(3.11) 



(3.12) 



3.1.1 Even metric perturbations 

Let us look for solutions of the tensorial wave equation fl3.8p . metric pertur- 
bations expanded in tensor harmonics according to [39] . Those of even parity 
can be written in the form 



^00 
hoi 

ho2 

hos 

/ill 

hu 
his 

h22 

h23 
h33 



dY 

He'' 



ime-'^^'-""''^hoir)Y{e), 
ye-*('^*-"*)if2(r)F(^), 

e-^^-^--^)hi{r)^, 
ime-'^'^'-""''^hi{r)Y{e), 



^2 ^-i{ujt-m<j)) 



K{r)Y{e) + G{r)^ 



dY 

'de 



cot eY{e) 



+G{r) 



-m^Y{9) + sin d cos 9 



dY 



(3.13) 



and metric perturbations of odd parity are instead found to vanish identically. 
Here, Y{6) is the 6'-dependent part of the spherical harmonics 11^(6', 0), 



40 



CHAPTER 3. THE TENSOR WAVE EQUATION 



solution of the equation fl2.10l) . It is convenient to introduce the notation 

K{r) = W,{r), Ho{r) = W2{r), H,{r) = Ws{r), H2{r) = W,{r), 
G{r) = W,{r), ho{r) = H~'We{r), h{r) = H''Wr{r), (3.14) 

in place of standard Regge-Wheeler-Zerilli notation for metric perturbation 
quantities. The wave equation (13. 8p . using the relations (I3.13p . (I3.14p and 
(I2.10p . leads to the following system of coupled differential equations: 



r'f f 



(2 - 3/)^ - ^[a;V + 2(1 - /)(1 - 4/) - fLjW^ 



2H\ dWi AH^r dWi 2H\L dW^ 
_^ 1 i _^ 2. 

f dr f dr f dr 

d^Wi 2(2 - 3/) dWA 1 



dr"^ rf dr r'^f'^ 

2H'^r dW2 4:H\ dWi 2rLH^ dW. 



[uj\^ - L/ - 10 + 6/ + 12(1 - /)2]iy4 



H ; ^ + 



f dr f dr f dr 



'2, 



r^/ r^f H r^f f 

^ = 7f 2/)^ - ;57jl- + 2/(3/ -2)11^5-;^, 

+ 4^[2 - (1 + /)L]iy5 - ^^\t^^ ^4 + ^W^i + ^1^2.(3.15) 
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To these equations one should add the de Bonder gauge components 



dWs iuL^^^ L iu^^^ f,,r 2(1-2/)^^^ 



dr 2/ ' Hr'^f ° 2/ " f \ " 2 J rf 

dW, 2(1-3/) ^^^ dW2 , dWi dW, 2L 

— — — VV4 H/— 5 L— W5 

dr rj dr dr dr r 

2L 4 2H^r 2iuj 

nr'^ r J j 

. !(^M,,^MH'.->.iiay.-;y.). (3.16) 

dr rf 2/ P 2/ 



Solving these two systems of coupled differential equations is clearly a 
hard task. A first progress can be done by introducing the same dimensionless 
radial variable seen in Chapter 2, that is 

X = Hr, (3.17) 

as well ClS db dimensionless frequency parameter 

n = H-^u. (3.18) 

A second step consists in reducing the set of second order differential 
equations to the normal form, i.e. considering the following rescalings: 

{W,,W2,Ws,W,,W,) = — =L=(/i,/2,^/3,/4,/5), 

xyl — 

{We,Wr) = (^tf,,^Jl-^y (3.19) 
It turns out that the system (13.151) can be cast in the form 



dx"^ 



Ai'fj, (3.20) 



where first order derivatives have been eliminated and with the only non- 
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vanishing coefficients listed below 

i_ x^n'^ + + Lx^ - 2 - L ^ 2 _ - 
Ai — 777- 7777; , ^1 — 77, Ai 



A, 



5 



A 



3 



A/ 



A,' 



— x^)^ ' 1 — x^' x^(l— x^) 

2L 

x^ 

,1 4 ,2 2x^ + x^L - L + x^n^ _ 3a;2 

A2 = , A2 = , A2 = , 

1 — x^' x^(l — x^)^ ' (1 — x^)^' 

A,' = —^, A,'- 



(1 — x'^y ' 1 — x^ 

, . 2^]x ,0 -x^ + ^]2a,2 ^ 3,2^ _ ^ _ 2 2^]x 

= A'i = A'i = 

^ (1-X2)2' ^ X2(l-x2)2 ' ^ (1-X2)2^ 

2L 



X2(l -x2)l/2- 



X2(l-X2)' (1-X2)2' (1-X2)2' 

, 4 -2x^ + x2L + x2^]2 + 5a;2 _ 4 _ ^ 2L 
A/ = — — , A 5 



x2(l — x2)2 ' ^ x2(l — x2)' 

4L 



X2(l -x2)l/2- 



4X^ + x2L - 5x2 ^ ^2^2 + 2 - L 



^ 5 _ -tu, -r ^ ^ uu. -r ^ -r ^ ^ A = - 

x2(l — x2)2 ' x2(l— x2)^/2 

,0 2 . 4x^ - L + x2^]2 + a;2^ _ 4^2 

Ai? = — — — -, - 



X2(1-X2)V2' ° X2(l-x2)^ 

2xn 



(1 -X2)2- 

A/ = ^ A.' = ^ A.' = 

^ x2(l-x2)l/2' 7 a;2(l _ 3,2)1/2' 7 x2(l-x2)l/2' 

2xr] 7_ -L + x2L + 4x2 + x2^]2 _ 4 

^7 - (l-x2)2' - X2(l-X2)2 • ^^-^^^ 



The rescaling (13.191) also implies for the gauge equation the following 
form: 



dx 

dx L dx L dx L dx^ 

£ = (3.22) 
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with the only non-vanishing coefficients hsted below 

R 1 ^ R 2 ^ R 3 — 1 

-D.3 = —Tz -Da = —TTTz -DS ~ 



(1-X2)' ^ 2(1 -X2)' ' X(l-X2)^ 

r4_ ^ p 5 _ R« ^ 



2(1 -a;2)' " 2(1 -x2)' ^ x(l-x2)i/2- 

2 r) 1 Q 



r1 r2 r3 

-D.c; — r-, -D.ts — 7Z -D, 



Lx(l-x2)' Lx(l-x2)' L(l-x2 

-D.c; — — -— r-, — r-, n> 



xL(l — a;2)' x(l — a;2)' x(l — s2)i/2 

^ 2x(l-x2)V2' 7 2x(l-x2)i/2' 7 2x(l-a;2)i/2' 
i?7^ = 7T^, 5r^ = --. (3.23) 

(1 — X^j X 



By differentiating Eqs. ( I3.22p and using into the new system Eqs. (EI 
as well as Eqs. (13. 22 p . gives a new system of ffist order equations 

ax ax 

^ = B,^ f, + C,^ + DJ„ 
ax ax 

^ = B,^f, + c/^ + D,U. (3.24) 

We list here the non-vanishing coefficients of this system. It is convenient 
to introduce the quantity 

A = 2x'^L - 2x^ + 2x2^]2 + 2 - L . (3.25) 

We have then 

B-^ = — (4x^L- 2 + x^n^ + 3x^-2x2^2 -x^-fiVL + L-Sx^L), 

x(l — x^)A 

B,^ = .^^{2 + 2x'Q' + 2x'L-L + 6x'-5x'), 
(1 — x^jA 

B2^ = — — ^— -(-4x2^ + 5x^^2 + 4x^L-2-9x^ + llx^ + L), 
x(l — x2)yi 

= -- — ^4TTT(-2x2fi2 + nV -3-3x^ + 6x2 + 3x^n2N 

(1 — x^yA 

2r]Lx2(l + fi2) 

Bo = 



(1-X2)3/2A ' 

, _ 6Lx(l-x2)(l-2x2) 

^2 - (1 _ ^2)3/2^ • (3-26) 
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4 ' 



(1-X2)A' 



= . — (-3x2- 2L-x^ + a;^Q2 ^a;^L + x2L + 2x^0^ + 4), 

a:(l — X'^jA 

x{i — x^yA 



= — ^^7tt7tt(-^^ + 2L + 2a;2L2 - Sx'^L + 20 VL + 12a;U - 12x'^ + 12x2), 

La;(l — x^)'-'^A 

B.' - ^^^3^(2L + .2^-4 + 4.2), 

^e' = -7^ ^^-^(-x2L-2x2 + 3x^L-L + 2 + QVL), 

(1 — x^yi'^A 

_ _^^(_2, + 2 + 2x2L + 2x202_22;2_202 + 02^), 
(1 — x^)^ 

= -—^^——{-L + l2x^ + 2-Ux^ + 2x^n'^ + 2x^L). (3.28) 
(1 — x^)A 



C2 



(2x2-l)(L-2) 

2(L-2)Ox(l-x2)V2 
LA 



(L-2) 



(3.29) 
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D2 = 7 ^^,r—(2x^-2x^ + 2n^-4:X^L + 7x'^L-6x'^n^ + n^x^L 

(1 — x^yA 

-n^L + Qx^n"^ + 2fiV - 3L), 
D4 = 7 ^-^r^(-i + 2x^ + 8x^ -6x'^ -x^L-2x^Q^ + 3x'^L-2x^Q'^ 

x{i — x^yA 

+ Qx^n'^ + 21] V - Sfi^x^L - 2x^L^ + + Sx^il^L + x^L^ - 2x^L), 

L(l — x-^yi-^A 

-L^ -8x^ + 8x^). (3.30) 

The final set of first order equations is then the union of Eqs. fl3.22p and 
(13.241) . implying for Eqs. (I3.22p the final form 



dx 

ax ax 
f = B,'f„ (3.31) 



where 



B, 



3 



4 



3x2 + L - 2 
xA ' 



LA ' 

x{5L + 4fi2) 
LA 



B5^ = — (-x2L + 4fiV-2 + x^L + 6x2-fiVL-4x^ + L), 

x[l — x'^jA 

B^^ = -— ^-^r-^—(-4nx^L + 4Lnx-2L'^nx + 2L^nx^), 

Lx[l — x^yi^A 

B^o = -T. ^^tt7^(6x' + 2x21]2 + 2x2L-8x2-L + 2), (3.32) 

x{\ — x^yi-^A 

C. ^ (3.33) 

and 

= ^, ^ r2 V?x^L + x^L^ - 8x^n'^ + 2L- 6x^L + 2x^L - L^). 

xL[l — x^)A 

(3.34) 



46 



CHAPTER 3. THE TENSOR WAVE EQUATION 




Figure 3.1: Numerical solutions for fi{x) with /i(O.l) = 1, f[{0.1) = 10, 
I = 1, ... ,7 and with i7 = 5 (left figure) and = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 



Such a first-order set of equations, fl3.24p and (I3.3ip . with a second-order 
equation for /i, that is 

= ^^A [L + ^-i^' + L + 3)x^] A - 2Ax' - 1)/, 

+ 2{x^-l)U-2L{x'-l)^f,}, (3.35) 

once used to replace derivatives into the system (13.201) , implies that the latter 
system is identically satisfied. 

At this stage, a numerical analysis of solutions can be performed easily, 
using the set of equations (I3.24p . (l3.3ip and (13.350 . 



3.2 Plot of the solutions 

We plot the solutions fi{x) — fj^x), beginning with /i(O.l) = 1, where i = 
1,...,7, and /((O.l) = 10. Figures from (3.1) to (3.7) describe the solutions 
for various values of I and Q. Figures from (3.8) to (3.14) describe the 
solutions for various values of I and fl, when /i(O.l) = 1 with i = 1, 7, and 
/{(O.l) = 100. 

As far as we can see, all solutions blow up (as the electromagnetic waves 
in Chapter 2) at the event horizon, corresponding to x = 1, since there are 
no static solutions of the wave equation which are regular inside and on the 
event horizon other than the constant one [351. 
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Figure 3.2: Numerical solutions for f2{x) with /2(0.1) = 1, f[{0.1) = 10, 
I = 1, ... ,7 and with ^2 = 5 (left figure) and ^2 = 10 (right figure). Increasing 
values of I correspond to more peaked curves on the right part of the plots. 




Figure 3.3: Numerical solutions for /^{x) with /3(0.1) = 1, f[{0.1) = 10, 
/ = 1, . . . , 7 and with f2 = 5 (left figure) and = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 




Figure 3.4: Numerical solutions for f4^{x) with /4(0.1) = 1, f[{0.1) = 10, 
/ = 1, . . . , 7 and with Q = 5 (left figTirc) and Q = 10 (right figure). Increasing 
values of I correspond to more peaked curves on the right part of the plots. 
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Figure 3.5: Numerical solutions for f^^x) with /5(0.1) = 1, f[{0.1) = 10, 
I = 1, ... ,7 and with fl = 5 (left figure) and fl = 10 (right figure). Increasing 
values of I correspond to more peaked curves on the right part of the plots. 




Figure 3.6: Numerical solutions for fQ^x) with /g(O.l) = 1, f[{0.1) = 10, 
I = 1, ... ,7 and with 1] = 5 (left figure) and = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 



/ \ 

/ 












Figure 3.7: Numerical solutions for fj{x) with /7(0.1) = 1, f[{0.1) = 10, 
I = 1, ... ,7 and with 1] = 5 (left figure) and O = 10 (right figure). Increasing 
values of I correspond to more peaked curves on the right part of the plots. 
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Figure 3.8: Numerical solutions for fi{x) with /i(O.l) = 1, f[{0.1) = 100, 
I = 1, ... ,7 and with fl = 5 (left figure) and ^2 = 10 (right figure). Increasing 
values of I correspond to more peaked curves on the right part of the plots. 





Figure 3.9: Numerical solutions for f2{x) with /2(0.1) = 1, f[{0.1) = 100, 
I = 1, ... ,7 and with f2 = 5 (left figure) and Q = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 




Figure 3.10: Numerical solutions for /^{x) with /3(0.1) = 1, /{(O.l) = 100, 
I = 1, ... ,7 and with Q = 5 (left figure) and Q = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 
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Figure 3.11: Numerical solutions for /^{x) with /4(0.1) = 1, f[{0.1) = 100, 
I = 1, ... ,7 and with fl = 5 (left figure) and Q = 10 (right figure). Increasing 
values of I correspond to more peaked curves on the right part of the plots. 





Figure 3.12: Numerical solutions for f^ix) with /5(0.1) = 1, f[{0.1) = 100, 
I = 1, ... ,7 and with ^2 = 5 (left figure) and = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 





Figure 3.13: Numerical solutions for JqIx) with /6(0.1) = 1, f[{0.1) = 100, 
I = 1, ... ,7 and with 1] = 5 (left figure) and O = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 
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Figure 3.14: Numerical solutions for fj{x) with ^(O.l) = 1, f[{0.1) = 100, 
I = 1, ... ,7 and with f2 = 5 (left figure) and = 10 (right figure). Increasing 
values of / correspond to more peaked curves on the right part of the plots. 
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Chapter 4 



Exact solutions in a de Sitter 
background 



4.1 Metric perturbations using the Regge- 
Wheeler gauge 

In Chapter 3, we have used the Lichnerowicz operator to derive the cou- 
pled system of differential equations for metric perturbations. However, it is 
possible to proceed in a different way. One considers the vacuum Einstein 
equations 

Rab - ^Rgab + ^9ab = 0. (4.1) 

If one introduces Qab = lah + e/^aft, where e is a parameter which controls 
the perturbation, one has 

Gl, + eGl, + A(7„5 + ehab) = 0, (4.2) 

where is the Einstein tensor with respect to the background 7afo and 
takes the form 

G\, = R\,-\{^abR' + K,R'), (4.3) 

where is the scalar curvature with respect to the metric 7a5, whereas 
and arc the Ricci tensor and the scalar curvature valuated to first-order 
in the metric perturbation hah- Thus, one has 

Gl, + A7a6 + e{Gl, + Khab) = (4.4) 
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and, since the background metric 'jab satisfies the vacuum Einstein equations, 
we must solve only the following equations: 

Kb - I ilabR' + KbR") + Ahab = 0. (4.5) 

4.1.1 Regge- Wheeler gauge 

Solving this system of coupled differential equations is clearly a hard task. 
However, we can introduce the Regge- Wheeler gauge |10] . With the notation 
used in (I3.13p . one has 



K{r) = W^ir), G{r) = 0, ho{r) = 0, hi{r) = 0. (4.6) 

Thus, on using the relations ^AO\f . ( KWf . ( KT7\f . f lSTTS]) and (gl]), the 
wave equation fl4.5p leads to the following system of differential equations: 



dWijx) 
dx 



dW2{x) 

dx 
dWi{x) 

dx 

WAx) 



+ 



-2x^ + (L - 4) x2 + 2 - L 

2x3 (a;2 - 1) 
{2- L-2VL^)x'^ - 2 + L 



WAx) + 



x'^{x'^ — 1) 



Woix) 







x^ — 1 
Wi{x) 

X 

W;ix) 
2-L 



2x%x^ - 11 
WAx) - 



WJx) 



x^ 



?^W,(x) + 4^WJx) 



x^ 



iL , , WAx) 

+ ^rT^W2{x) + , r 



2nx'^ 



X[X^ 



W,{x) + '-^^^^^W.ix) 



+ 



2x3 
(L + 2r]2)a;2 



2^]x2 

2-L. 



2x3 (x2 



1 



-Wi[x) 



(4.7) 
(4.8) 

(4.9) 
(4.10) 

(4.11) 



At this stage, we can solve this system exactly in terms of the Heun 
general functions (see [21] and Appendix D). In fact, if we use the equation 
(14. lip , we can isolate Wi{x) and, by using the equations (14. 8 p and (14. 9p . we 
obtain a second-order differential equation for W^i^x), that is 



WAx) = 



(4.12) 
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where A{x) and B{x) take the form 



4(L + 2 Q^)x^ + 2L{L - 2)x'^ - L{L - 2) 
^"^^ ~ xix^ - 1)[2(L + 2 Q2)^2 + _ 2)] 

^ 2 r)^(I + 2 r)^).r^ + 1(1 - 2)(I + r)^).r^ -LHL-2) 
^' x2(x2-l)2[2(L + 2 02)x2 + L(L-2)] ^ 

Eventually, we have 



W4{x) = (x^ - 1)2" [Cix ^ ^/feunG'(ai,5i,Q;i,^i,7i,5i, x^) 

+ C2x'-HeunG{a2,q2,a2,f32n2,S2,x'^)] (4.14) 

where 

/(2 + Z)(l-/2) 



?1 



2(2Q2 + /2 + 

(/ - + 2){i + - tni + 4ii] + 3/ + 2/^ - 1^) 

8{2n^ + P + 1) 



A = ^(iQ-Z-1), 

5i = + (4.15) 

and 

/(2 + 0(l-/') 



02 



?2 = 



2(2Q2 + /2 + ;)' 

^(1 - 0(1 + / + + 3il]/ + Qin + 41 + 5P + /3) 

8{2Si^VPVf) ' 



a2 = ^(iQ + Z), 

P2 = ^(^^) + / + l), 

3 , 

72 = 2+^' 

^2 = iQ + 1. (4.16) 
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At this stage, if we use the solution for W^i^x) and the equations fl4.8l) and 
(14.91) . we obtain the solution for W2{x). Moreover, by using the equations 
fl4.10p and (14. lip , we have the solutions for Wi{x) and W^^x). 

The decoupling technique and analytic formulae are original [21] but in 
this context, we want to solve explicitly the equation (14. 5 p also in the de 
Bonder gauge (see Chapter 3) and, eventually, to compare this result with 
that obtained in this Chapter by using the Regge- Wheeler gauge. 

Moreover, other powerful techniques are available for solving a tensor 
wave equation in curved space-time. The interested reader is referred, for 
example, to Appendix E. 



Conclusions 



Gravitational waves are been here considered as metric perturbations in a 
curved background rather than in the flat Minkowski one, since several in- 
teresting physical situations can be discussed in this framework. 

In Chapter 1, we have seen that in the case in which the de Bonder gauge 
is imposed, its preservation under infinitesimal space-time diffeomorphisms 
is guaranteed if and only if the associated covector is ruled by a second-order 
hyperbolic operator. Moreover, since in this case the Ricci term of the wave 
equation has opposite sign with respect to the wave equation of Maxwell 
theory, in the Lorenz gauge, it is possible to relate the solutions of the two 
problems |25j . 

In Chapter 2, we have completely succeeded in solving the homogeneous 
vector wave equation of Maxwell theory in the Lorenz gauge when a de Sit- 
ter space-time is considered. One component of the vector field is expressed, 
in its radial part, through the solution of a fourth-order ordinary differen- 
tial equation obeying given initial conditions. The other components of the 
vector field are then found by acting with lower-order differential operators 
on the solution of the fourth-order equation. The whole four-vector poten- 
tial is eventually expressed through hypergeometric functions and spherical 
harmonics. The decoupling technique, analytic formulae and plots are com- 
pletely original [20] . 

In Chapter 3, we have extended this method to the wave equation of met- 
ric perturbations in de Sitter space-time and we have written the Lichnerow- 
icz operator P^^^'^ which results from the expansion of the action functional 
to quadratic order in the metric perturbations. The basic equations of the 
theory read therefore as 

Pab^'h^ = 
^a{h) = 
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where $0 is the supplementary condition for gravitational waves. Eventually, 
a numerical analysis of solutions has been performed. 

In Chapter 4, we have solved explicitly the Einstein equations for metric 
perturbations on a de Sitter background. In fact, by using the Regge- Wheeler 
gauge, the coupled system of differential equations to first-order in the metric 
perturbation hab is been solved in terms of the Heun general functions. 

In future work, we want to solve explicitly the equation (14. 5 p also in the 
de Bonder gauge (see Chapter 3) and, eventually, to compare this result with 
that obtained in Chapter 4 by using the Regge- Wheeler gauge [2T] . 
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Appendix A 
Bivectors and biscalars 



In Eq. fll.SSp . g°'^, is the geodesic parallel displacement bivector (in general, 
bitensors behave as a tensor both at x and at x') which effects parallel dis- 
placement of vectors along the geodesic from x' to x. In general, it is defined 
by the differential equations 

a^" = = 0, (A.l) 

jointly with the coincidence limit 

\m,g^,=\g^J\ = (A.2) 

x'— >a; L 

The bivector when acting on a vector B'^' at z', gives therefore the vector 
^ which is obtained by parallel transport of B'^' to x along the geodesic 
connecting x and x', i.e. 

B" = g"^, B''. (A.3) 

In Eq. fll.56p . A(x,x') is a biscalar built from the Van Vleck-Morette 
determinant 

D{x,x')=det{a,ab') (A.4) 

according to 

A (x, x') = ^ D{x, x') ^ (A.5) 
V-7(a;) V-7(a;') 

The biscalar A{x,x') has unit coincidence limit: [A] = 1; as a function 
of X (resp. x'), it becomes infinite on any caustic formed by the geodesies 
emanating from x' (resp. x). When A diverges in this way, x and x' are said 
to be conjugate points [H]. 
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Appendix B 
Derivatives of /q 



The higher-order derivatives of /o in Chapter 2 get increasingly cumbersome, 
but for completeness we write hereafter the result for fdix), i.e. 



fo{x) 



4ai(ai + 1)61(61 + 1) 



a;'+'(l - a;2)"*^F(ai + 2, 61 + 2; di + 2; x^) 



+ x^-\l-x'') 



di{di + l) 

^2aibi 



di 



X 



X [(2/ + 1)(1 - x^) + 2inx^]F{ai + 1, 61 + 1; di + 1; x^ 



+ x'-\l-x')- 



\n-2 



l{l - l){x^ - If - ^(^2 - l){lx + 2{l + l)x') 



+ (2iO-02)x^ 



+ i C'a ^ (^4, ai ^ ai + 1, 61 61 - 1 



(B.l) 
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Appendix C 
Lichnerowicz operator 



The Lichnerowicz operator on metric perturbations is obtained by expansion 
of the Einstein-Hilbert action to quadratic order, i.e. 



6^S 



EH 



where 



S 



1 



EH 



Sgf"" = 

-4 



IGtc 



d x^/g R 



M 



and g = —det{g^y). The first functional derivative of the action is 



5S 



d'x^-^R+ [ d^x^^4^RP''+ [ d'x^g,, 
Sg'"' J ^ 6g>''' J ^ ^ 



SRP" 



where the first two integrals can be written as 



d^x^-/lR+ fd^xVg^'^ 



6g 



d^Xy/g 



(C.l) 
(C.2) 

(C.3) 

(C.4) 



Thus, the first functional derivative provides the vacuum Einstein equa- 
tions. In fact, the integrand of the third integral of (1C.3P is a total derivative, 
that is 

gp<rj^ = -V,[5^V.5(x,x') + 6y,6{x,x') - g^y<'6{x,x')\. (C.5) 

At this stage, it's possible to write the second functional derivative of the 
Einstein-Hilbert action, i.e. 



bg^'bg 



Sg^'Sg 



d^x ^ X R+ / d^'x-^^R'"' + / d^'x 



Sgi"' Sg 



■ye 



Sy/g SRT 
5gP-^''' 6g 



7£ 



+ / d^X 



5gie §gfMU 



Rp"" + [ dS^ f^'/ RP"" + [ d^x^ 



5g,a SRp'^ 
Sgp" 5g^^ 



(C.6) 
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and if we note that 



= ^y/g 9je9nu - ^V9i9nj9ue + 9ne9v'^), (C.7) 



one has the expression for the first integral, i.e. 

1 

Sg-yeggHu ^ = lV9i9fil9i^e + 9ne9uy - 9ie9nu)R, (C.8) 

moreover, one has 

l^l?^"'^ = -lV^9,A9n9ae + 9,e9a,)R'"' = -\^99,.R,e (C.9) 

and similarly, the fourth integral is 

J^eJ^uR'" = -\V99,e{9p,9a. + 9pu9.,)R'^ = -\V9 9,eR,u. (C.IO) 

At this stage, we note that the variation of the Ricci tensor respect to 
the metric can be written as 



1_„_ 1 

4 



-{5f;5l + <)T5^)V,V" - -^..VV^ (C.ll) 



Thus, the third integral is 
5^ dRf"" _ 1^/1 



SR^"^ 1 f 1 1 

V'^^'"^ = -l^V9\l^9A'^l^e + VeV,)-9,v9,eVcN''^ (C.12) 

and for the fifth integral, if we note that 

^gie^gixv ~ 29(Tv{,9pl9tie 9pe9lJ.l) 29plJ.{.9<y^9i'e 9<T£9v'l) 

+ ]j^9ai,{gp^gve + gpegv-y) + ]j^9pu{9<T-i9p^ + gaegfi-f), (C.13) 

one has the expression 

V9 ^gJ^^^^ R^" = 2 [5/.£^i^7 + 9^'rReu + 9ueRJ^l + 9ujRt,e]- (C.14) 
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Moreover, for the last integral, one has 

V 9 ^^^^ ^^^j 2 V 5' [9fJ.p9au + gpu9crti) ^^^^ 

= ^V9 [9ep'^^'^u + Qy^i^e'^p + S'e^.V^V^ + g^fpV eV ^ 

- {9ip9eu + gepg^v)VaV'' - ^^.(V^V, + V,V^)]. (C.15) 

At this point, adding the terms obtained, one has the following operator: 

1 

'^pv-ie = -TV^ [{g-ypgeu + gepg^u " '^gpug-^e)^ cN" 



- g^p^j^iy - giy-y^e^p - - g^pV u 

~l~ 2((7^i/-R'y£ ~l~ g-ye^pu gpe^wy gp-y^eu g^e^-yp gu-i-^pe) 

- {gp^gue + gpegu-y - g^egpu)R]- ( C . 1 6 ) 

However, T^^-ye is a singular operator. One has a supplementary condition 
for the Sgfj^^, that is 

PrSgpu = o. (C.17) 

Thus, using the supermetric G^^^"^ , one has 

P'^"" = -^G^''''"(A)g^,°, (C.18) 

with 

Qp." = -(^pV. + (5^V,) (C.19) 

and 

G'M.p.(^) = ^{g^'g"" + g'^'g"" + Xg'^'gn. (C.20) 

where, however, not all the values of A G R are acceptable. In fact, the 
supermetric G^^'"'{X) must be invertible, that is 

G^^'"^{\) G,.er{f{X)) = \{5^rK + 5^5.^), (C.21) 
with ^ 

GpaerUW) = l^igpegar + gprg^e + /(A)5'pafl'r£] • (C.22) 

Thus, using KVl^ and (ESI into flCTzTj) . one has 



G"''^^'^(A) G,.e.(/(A)) = \{5^J: + 4^5^^) + ^^'^'^^7. 



A /(A) nA/(A) - 

2 2 4 

(C.23) 
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where n is the dimension of the space-time. Now, if we want that (1C.23P is 
equal to (IC.21|) . we must impose that 



A /(A) nXfjX) 
2 2 4 



Thus, solving this equation for /(A), one has 

2A 



(C.24) 



/(A) 



2 + nA 



(C.25) 



which is the coefficient of the inverse supermetric. Eventually, we note that 
the supermetric is invertible if 



A^— . 



n 



(C.26) 



At this stage, using (1C.19P and ( 1C.20|) . one has 



'-G>"'f"'{\)Q 



per 



and, eventually, the supplementary condition becomes 



(C.28) 



Now, we compute the supplementary condition applied on the metric 
perturbations, i.e. 

$°(V,A) = P'^^"(A)V, (C.29) 
where h^u = 5gfj,y. Thus, using (1C.28|) . becomes 
1, 



[g^'-V'^K. + g'^^V^K. + Xg^'^V'^KA 

1, 



V„ ( /i''" + -g^'^gp^h 



per 



(C.30) 



where if we put A = — 1, one has the de Bonder's gauge. 

At this stage, we can to compute the invertible operator S^^u^yei that is 



T =9 — P "'o »P ^ 



(C.31) 
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Thus, if we note that 



one has 



^M^" = ^('^^V. + + \g,.V-), (C.32) 



- ^AV^?7.V„V". (C.33) 
Now, using ( IC.16P and fIC.SSp . the operator J^^^e becomes 

J^lU'ye{>^) = -^Vg {[g-/f,geu + ge^lg■yu - (2 - A^)5(^^5(^£]VaV° 

+ (1 + A)(7^,(V^V, + V.V-,) + (1 + \)g^e{^^Nu + V,V^) 

- ^M7(V.Ve - V,V,) - g^,eiSuV^ - V^V,) 

- gu^i^tNe - V,V^) - g^ei^^V^ - V-,V^) 

~l~ '^i.g ^ii^Pje ~^ g^yePfiu gfiePvy g^-yPeiy gi/e-PjiJ, gvyPfie) 

g,eg,.)R}. (C.34) 

At this stage, we note the following relations: 

-^M7(V.V, - V,V,) - ^;..(V,V^ - V^V,) - ^.^(V^V, - V,V^) 

-fl'i/e(V^V^ - V^V^) = gr-y{RJ,^^ + Rl^J + g^r{.Rl^ue + Rlv-i) 
~^gT£{R^u-y ~^ Ryfi^f) ~^ QuriR-y^e ~^ Refir) ~ '^Rfi-yus + '^Rfieu'y- (C.35) 

Eventually, with A = — 1 and using fIC.SSp . one has 



1 



3^ fiuje = -^y/g{ [g^^^geu + get^g-fu - gfj,ug-(e] (V^V" - R) + 2R^^^^ + 2R^ 



+ '^{gi-iuR'ye + g-reRfiu ~ g^ieRwr " giM-yReu ~ gueR-y^i ~" gwyRfie)} {^.36) 
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Appendix D 
Heun functions 



The five multi-parameter Heun equations have been popping up with sur- 
prising frequency in apphcations during the last 15 years. Heun equations 
include as particular cases the Lame, Mathieu, spheroidal wave, hypergeo- 
metric, and with them most of the known equations of mathematical physics. 
Five Heun functions are defined as the solutions to each of these five Heun 
equations, computed as power series solutions around the origin satisfying 
prescribed initial conditions. 

The General Heun equation, with four regular singular points in {0, 1, a, oo}, 

is 



+ 



7 5 

- + 

z z — 



+ 



a 



d 



+ 



a/3z — q 



dz z{z — l){z — a) 



y{z) = 0. (D.l) 



The solution to this equation is implemented in this work as the HeunG 
function. The sum of the exponents of the singularities of Heun's equation 
is equal to two and the parameter e is expressed in terms of the other ones 

by 

e = Q; + /3 + l-7-(5. (D.2) 

The other four Heun equations are confluent cases, obtained from the 
general Heun equation above through confluence processes. The Heun Con- 
fluent equation, with regular singular points in {0, 1} and irregular in {oo}, 
is 

5 \ d q — a/3 q 

z 



dz^ 



+ 



z z — 1 



— e 



dz z — 1 



0, 



(D.3) 



having for solution the HeunC function. The Biconfiuent equation, with 
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regular singular point in {0} and irregular in {oo}, is 



/ ^ l + a\ d ^ 1(1 + a)/3 + S] ^ 

(D.4) 

having for solution the HeunB function. The Doubleconfluent equation, with 
two irregular singular points in { — 1, 1}, is 



dz^ 



a + 2z + az'^ - 2z^\ d 5 + {2a + ^)z + ^z'^ 



{z^ - ly 



dz ~'~ 



(^2 _ If 



viz) = 0, (D.5) 



having for solution the HeunD function. Eventually, the Triconfiuent equa- 
tion, with irregular singular point in {oo}, is 



y{z) = 0, 



(D.6) 



having for solution the HeunT function. The Heun functions, HeunG, He- 
unC, HeunB, HeunD and HeunT, arc defined as the solutions to the corre- 
sponding General, Confluent, Biconfluent, Doubleconfluent and Triconfluent 
Heun equations. These solutions are constructed as power series solutions 
around the origin, for certain initial conditions. 

The power series solutions at the base of the functions' definitions have 
restricted radius of convergence in the HeunG, HeunG and HeunD cases, 
where the numerical evaluation is done using analytic extensions, exploring 
closed form identities satisfied by these functions, as well as series expan- 
sions around different singularities. For arbitrary values of the parameters, 
however, closed form formulas for the connection constants relating series 
expansions around different singularities, are not known. 



D.l Some important facts about Heun func- 
tions 

• The coefficients entering the series expansions represented by the Heun 
functions satisfy three term recurrence relations. A solution to these 
recursion equations is not known in the general case, so a closed form 
for the series's coefficients is not available and the computation of - say 
- the nth coefficient requires the explicit computation of all the previous 
ones. 
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• They are more general than the rest of the functions of the mathe- 
matical language in that they contain most of them as particular cases. 
Consequently, the Heun equations cannot have their solution expressed 
(but as infinite sum power series) without using the corresponding Heun 
functions. 

• The Heun functions have a rich structure and so satisfy a rather large 
number of identities. 

• Because they have such a rich structure and include as particular so 
many functions, including the Mathieu, Lame, Spheroidal Wave and hy- 
pergeometric functions, the interrelations between them and the Heun 
ones are a source of many nontrivial identities between the former. 

• Due to the enlarged structure of singularities (if compared for instance 
with hypergeometric functions) the Heun functions are increasingly ap- 
pearing in the modeling of different types of problems in applied math- 
ematics. 
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Appendix E 



Wald's method on a de Sitter 
background 

Depending on the coordinate system used, there are many ways of viewing 
de Sitter space. It is has been studied largely because of the central role 
it plays in almost all inflationary scenarios of the early universe. Roughly 
speaking, the expansion is driven by a large cosmological constant which 
appears due to the energy density of a false vacuum. In this regard the 
description of de Sitter space by Robertson- Walker (RW) coordinates has 
tended to be the natural choice for most workers because of their obvious 
cosmological significance. In this coordinate system, constant time surfaces 
appear homogeneous and isotropic. 

However, yet another picture of de Sitter space time is provided by static 
coordinates. Here homogenity of constant time surfaces is lost but space- 
time appears static within a horizon distance, a very different state of affairs 
from that prevailing in the RW description. Gal'tsov and Nunez [42j have 
studied gravitational field perturbations in a de Sitter background described 
by static coordinates. To accomplish this its have employed the technique 
of Debye potentials introduced by Wald [43J . This treatment is particularly 
attractive because it reduces the problem of solving the sourceless equations 
for fields of different spin s to that of solving a single differential equation 
for the Debye potentials with free parameter s. 
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E.l Wald's technique for the Debye poten- 
tials 



Wald has described ^3] a remarkable technique for solving field equations. 
This technique requires writing the differential equations in two ways, which 
are related to obtain a third equation for the so-called Debye potential. The 
solution of the original field equations are expressed as operators acting upon 
this potential. 

Explicitly, Wald's method takes the field equations as its point of depar- 
ture 

\s\^ab^'' = 4:7ris\ja, (E.l) 

where s is the spin weight, \s\^ab are field operators, ip^ is the "middle" 
potential and \s\ja is the source for the field. Next Teukolsky's work [H], 
which makes use of the Newman-Penrose formalism [15] to express the field 
equations, is then incorporated and the field variables ip are expressed as 
operators acting upon the middle potentials 

siJ = sMaif". (E.2) 

Wald went on to obtain a second way of writing the field equations 

^,n,M„y;'^ = 47r,r'^|,|j„ (E.3) 
pp* 

where p is a spinor coefficient, is a Teukolsky operator, ^r" is a source 
projection operator and the other quantities are defined as in Eqs. ( 1E.1|) and 
(0. 

By operating on Eq. (lE.ip with ^r", it is possible to equate the projected 
left-hand side of Eq. dH]) with that of Eq. (iKSj) 

,r%|£,6/ = ^,n,M,/. (E.4) 
pp* 

From here Wald proceeds to the operator identity, 

sr%|£afe = ^.D.Mb. (E.5) 
pp* 



However, this step needs some justification. For example, in the case 
of the scalar function, V^(y9 = 9^<^ does not imply the operator identity 
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= dfj_. The validity of equation Eq. (lE.Sp has, however, been estabhshed 
for a wide class of space-time which includes the de Sitter space |16] . 

The adjoint tilde of the operator identity ( lE.Sp is now taken. By making 
use of the fact that the operator \s\^ab is self-adjoint for the fields of physical 
interest, we obtain 

If the function ^H, satisfying 

^ ^ ^ =0, (E.7) 



is introduced, then on account of Eq. (lE.6p . also satisfy the equation 

\s\^ab sT""" -s2 = 0, (E.8) 

sE is known as the Debye potential. For integer spin fields, the operator \s\^ab 
is real, in which case {gf"- _sS)* also satisfies Eq. flE.Sp . For the spin ^ field, 
however, this is not so. 

In the integer spin case, we define 

s^"^^^ = ^{sf'' -.2 ± i,f^ _,S)*}, (E.9) 
where a = 1 for (+) and a = i for (— ). So 

\s\8,ab s^"^^^ = 0, (E.IO) 

that is, the function sV''*'-^^ satisfy the sourceless fields equations flE.ip . For 
the spin ^ field, the solutions are given by f lE.Qp without taking the complex- 
conjugate. Thus, the problem of solving the field equations reduces simply 
to solving Eq. (]E.7p . for the Debye potential. 



E.2 Wald's method in de Sitter space 

In a cosmological context it is natural to describe de Sitter space in Robertson- 
Walker (flat) coordinates 



ds^ = de - e2"*(rfr2 + r^de^ + sin^ 9 dip^). 



(E.ll) 
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However, Wald's method is most directly applied in static coordinates 
(r, X; 6', if), where formally de Sitter can be treated as a member of the Kerr 
family. These are related by 

r = t- — ln(l -aW"*); x = ^e"*; 6 = 6; <f = ip. (E.12) 
2a 



In this coordinates the line element takes the form 

A, 2 X' 



ds^ = ^dT^- ^dx" - ide'' + sin^ edip^), (E. 13) 



where 

A = x'(l-aV) (E.14) 
and a is related to the cosmo logical constant A by = |-. The de Sitter 
horizon is given by x+ = o^ '^- 

We now introduce the Kinnersley null tetrad given by 

A 

A / y2 



P = 1^,1,0,0 



rrf = -^{0,0,l,icsc9) . (E.15) 

With respect to this choice of tetrad the nonzero spin coefficients are 

1 A A' , , 

^ = ^^2^^' ^^^^^ ^ 

and 

where the prime represents differentiation whit respect to x- We notice also 
that the Weyl tensor vanishes; thus in the Petrov classification, de Sitter 
space is a type N space (a particular case of type D space). 

We work with the operators Dp, D^, Lg and L+, which where introduced 
by Chandrasekhar |17] in the black hole context. In the case of de Sitter 
space, they are related to the Newman-Penrose operators D, A, 6 and 6*, in 
the following way 



A' Ao^ A' 

Dp = D+p-; -^D; = A-2p— (E.18) 
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and 



-^L+ = 5 + 2g/3; --^L, = 5* + 2g/3, (E.19) 



or, more explicitly, 



A' 

h — 

A A 

A A 

'C'q = dg + qcotO — icscOd^p, 

£+ = (E.20) 



Dp = d^+p— + —dr, 



With respect to the adjoint operation, 

= -V^, (E.21) 

so that the adjoint Chandrasekhar operators are 

Dp = -p^D^pp^, 
D; = -p'Dtpp\ 

^+ = -^1-,. (E.22) 

The derivation of the Debye's potential has been already done for a large 
class of type D and type N space-times, which include the de Sitter space as 
a particular case of type N space [38], so we give just the results for the field 
perturbations in the absence of sources in terms of the corresponding Debye 
potential 

For the gravitational case s = ±2, the field perturbations ±2^^^ are given 

by 

±2/^^^^^) = ^{±2f''%2^ ± c.c.}, (E.23) 



where 



+2f/^- = FF^J:^/:2 + l^''m*''^V2p-^/:2DoP^ 

4rn .3 



+ mV^D+p-^2)+pl (E.24) 
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For the electromagnetic case s = ±1, the vector potential ±iA'^ is given 

by 

±i^'^^^^ = ^{±ir^iS±c.c.}, (E.25) 

where 

^if" = n^f^Lt -m^'^^V+p. (E.26) 
V 2 2 

The Debye potentials for each field satisfy a second-order differential equa- 
tion, which leaving the spin projection s as a parameter can be written as a 
single master equation, obtained from the Teukolsky equation 

.□%S = 0, (E.27) 

on account of the identity 

sB? (E.28) 
where the Teukolsky operator, is given, for s > 0, by 

, □ = ADiZ}+ + - 2{2s - l)xdr - 2(2s - - l)aV (E.29) 
and, for s < 0, by 

, □ = ADt^.Do + Li+s^t, -2{2s + l)xdr - 2{2s + l)a V- (E.30) 

The explicit form of the operators ^r" and sMa, needed in these derivations, 
are given in |12]- Because we are working in static coordinates, we expect a 
time behavior e~*'^^; we observe, moreover, that the Debye potentials factor- 
ize into radial, angular, and temporal parts as follows: 

= j sR{x) -sYie, ^)e-'^^du. (E.31) 

Furthermore, it is clear that equation f lE.27p separate into radial and an- 
gular parts. Substituting (iKSTj) into IK27} and by using (E:29|) and f EIMj) . 
we obtain for the radial functions {s > 0) 

[ADiDj + 2{2s - l)tujx - 2{2s - l){s - l)aY] sR{x) = sXsR{x) (E.32) 
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and, for s < 0, we have 

[AD+,,Do + 2{2s + l)tujx - 2{2s + l){s + 1)«V] sR{x) = sXsR{x)- 

(E.33) 

For the angular functions, we obtain 

where <jA is some separation constant. Eventually, it is possible to show that 
for the radial functions, one has |12] 

sRix) = sCiX'''-'f(i + s + 1,-1 + s,1 + s+'^^^'^^ 



2a 2ax 



X 



l+_ax 
1-ax 



2q;' 2ax 

_ iui 
\ 2a 



X |1±^) . (E.35) 
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